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Q ■ Abstract. N. Kishore, Proc. Amer. Math. Soc. 14 527 (1963), considered 

the Rayleigh functions an{i^) = Y,^=iJvk'^i n = 1, 2, . . ., where rLj^k are the 

(Sj \ (non-zero) zeros of the Bessel function Ju{z) and provided a convolution 

type sum formula for finding a^ in terms of o"i, . . . , (Jn-i- His main tool was 

^ ■ the recurrence relation for Bessel functions. Here we extend this result to a 

^ ■ larger class of functions by using Riccati differential equations. We get new 

r^ . results for the zeros of certain combinations of Bessel functions and their 

'j^ \ first and second derivatives as well as recovering some results of Buchholz 

for zeros of confluent hypergeometric functions. 



^ • 1 Introduction 



The Rayleigh functions are defined, e.g., in |I|, p. 502], by the formula 



oo 



where ±juk are the zeros of the Bessel function 



e:: a„,(z.) = 5:j-2", n=l,2,..., (1) 

s 

5^ \ They form the basis of an old method due to Euler, Rayleigh and others for evaluating 

the zeros. For example, in case z/ > — 1, the inequalities 

WnivT^'"" < jIi < a^{y)lon^^{y), n = 1, 2, . . . 

provide infinite sequences of successively improving upper and lower bounds for j^^. 
Several authors have considered the question of finding "sum rules" or formulas for 
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an{i^)- By a method originating with Euler (see |I|, p. 500, ff.] for details; various 
ramifications were considered recently in p[), we can find all the (Jn{v) in terms of the 
coefficients in the series @). If we want to deal (as in |^) with properties of the o"„(z/) 
as functions of z/, there is a useful compact convolution formula due to Kishore [ffl 



1 n— 1 

(ynii') = Yj (yk{i^)(yn-k{y), (3) 

^ ^ "^ k=l 
from which the 0"„(z/) may be found successively, starting from 

ai = l/[4(z.+ l)]. (4) 

The question arises whether there are Kishore-type formulas for sums of zeros of other 
special functions such as the first and second derivatives of the Bessel function. In [^ 
there is a variant of this result for the zeros of the more general function 

N^{z) = az^j'liz) + hzj'^{z) + cJ^{z) (5) 

considered by Mercer 0. The result of gave a method of finding the reciprocal 
power sums 

oo 

rn{y) = Y.^~ul''^ n = l,2,.... (6) 

fc=i 

where x^k are the zeros of the function N^{z). The main result of expressed r„ in 
terms of r^, A; = 1, . . . , n — 1 and ak, k = 1, . . . ,n. It seems desirable to express r„ in 
terms of r^, k = 1, . . . , n—1 only. We do this here by using the Riccati equation satisfied 
by z~^/^A^^(x^/^). We record also the second order linear differential equations satisfied 
by N^(z) and by z~'^/'^N^{z^^'^) since these do not seem to appear in the literature and 
may prove useful for other purposes. 

In §4, we apply the same method to get power sums for zeros of confluent hypergeo- 
metric functions. 



2 Differential equations for functions related to 
Bessel functions 

The Bessel function y = Jy{z) satisfies the differential equation 

zy + zy' + {z^ - u^)y = 0. (7) 

and the function y = zJ^{z) + cJu^z) satisfies [0, p. 13] the differential equation 



Here we record the more general second order linear differential equation satisfied by 
the function 

y = N,{z) = az^jl{z) + hzj'^{z) + cJ,{z). (8) 



It is 
where 



z'^y" + A{z)zy' + [B{z) + z^ - u^]y = 0, (9) 

— 3a^z^ + pz"^ + q 



A{z) 
Biz) 



2a{a + b)z^ + 2rz^ 
a?z'^ — pz"^ + q 



with 



and 



p = 2a{ai/ + c) + (a^ - 6^), 
q= (ai/^ + c)^ -v'^{a-hf, 



r = az/^(3a — b) + c{a + b). 
We found the equation (|^) by repeated use of 

zJl{z)=uJ,{z)-zJ,+,{z) (10) 

to express the derivatives J'^\z)^ n = 1, 2 ... in terms of Ju{z)^ Jy+i{z) and discovered 
an appropriate vanishing linear combination of N^{z), Nl{z), N"{z). Of course, once 
(H) is known, it is easy to verify that N^{z), given by (||), satisfies it. 

It is convenient to consider the function 

y,{z) = z-^I-'NM''') (11) 

where we choose that branch of 2^^^ which is positive for z > 0. Using (^, we find 
that the function y^iz) satisfies 

4t2^ + [4z/ + 2 + 2A(i'l'')\t^^ + [t - z/ + vAit^l'') + B{t^l'')\y = 0, (12) 

It is well known that if y satisfies 

y" + P(t)y' + Q{t)y = 0, (13) 

then u = y'/y satisfies the Riccati equation 

nil 

— + P{t)u + Q{t)+u^ = 0. (14) 

Applying this to ([T2|), we find that, with y^iz) given by (H), u = yl{z)/y^{z) satisfies 

du 



At{aV -pt + q) 



^ + " 



+ 4[a^(z/ - l)t^ - upt + q{u + l)]u + 
+ aH"^ + b + 4aV - 2a{a + b)]t + 2i^p + q + 2r = 0. (15) 



3 Functions of Rayleigh type 



The even entire function z ^N^i^z) has an infinite set of zeros ±t„, n = 1, 2, . . . with 

so the zeros of yv{z) are Cfc = ^i, with 

ElC'l<oo. 
Thus 



y,{z) = z'^/'N^iz'/') 



au + c+ {b — a)iy 



fc=i V Cfc/ 



(16) 



The constant muhiphcative factor is got from the series (0). The vahdity of this infinite 
product expansion follows from facts on entire functions of finite order |^, Ch. 8] . 

We may differentiate (|T6D logarithmically [|], to get 

yliz) ~ 1/a ^ 1 oo .n 



Vuiz) 



Sle 



k=l 



z/Ck 



Y-Y—- 

k=l Sk n=0 Sk 



This gives 



2z^ 



yu[z) 



-2 E E ^VC- 

k=l ra=l 



But we may interchange the orders of summation here (since the iterated series con- 
verges absolutely) to get 



2z 



yu{z) 



oo oo 



-2 E ^" E C" = -2 E ^nZ"", 
n=l fc=l n=l 



(17) 



where 



Using 



we get 



oo 
^n = E C"- 



U 



U 



E 

k=2 



E ^fc+1^ 



fc=0 



fc-l 



/ , T~mT~k~ 



.m=l 



^k-2 



Substituting in ([TsD , and comparing coefficients of powers of z, we get 

2z/p + q + 2r 



n 



4g(z/ + l) 



4g(z/ + 2)t2 = Aqr'} + Aupn -p- 4aV + 2a{a + b), (19) 

4g(i/ + 3)r3 = 4p(z/ + 1)t2 - ^c?{y - l)ri + a^ + SgriTs - 4pr2, (20) 

and, for A; > 3, 

g(A; + z/ + l)rfe+i = p(A; + z/ - l)rfc - a^(A; + z/ - 3)rfc_i 

k k-\ k-2 

~\~<1 2^ 'T'm'Tk-m.+l ~ P 2^ 'Tm^k-m + O- /^ 'T-m'Tk-m.-l (21) 
m,=l m=l m=l 

In the special case a = 6 = 0, c = l (and hence p = 0, g = 1, r = 0), where we are 
deahng with the zeros of the Bessel function, these reduce, as they should, to (^ and 
the convolution formula (0) for cr„, n = 2, 3, . . .. 

In the special case a = c = 0, 6=1 (and hence p = — 1, g = — z/^, r = 0), we are dealing 
with the non-trivial zeros of the function Jl{z); (|^), (^^ and ( pi] ) become 

4(z/ + l)z/ 

_ -4z^Vf - 4z/ri + 1 
^'" -4z/2(z/ + 2) ' 

z^^(z/ + 3)r3 = (z/ + l)r2 + 2z/Vir2 - r^^ 

and for A; > 3, 

- z/^(A; + z/ + l)rfe+i = -{k + u - l)Tk 

fc-l k-2 

m=l m=l 

In particular, these lead to 

^r-/ 1-4 1 Z/2 + 8Z/ + 8 ,^^, 

^' = t}''^'^ =I^^V+l)V + 2)' ^^^^ 

= y^f.. 1-a ^ 1 Z.3 + 16Z.2 + 38Z. + 24 
3 Z^Uukl 32z/3(z/ + l)3(z/ + 2)(z/ + 3)' ^ ^ 

the same results as are obtained by the power series method in 0. 



4 Confluent Hypergeometric Functions 



Buchholz [|To| studied the nontrivial zeros a\ of the function 



zb/2^-z/2 



and showed that these zeros are all simple and that there are infinitely many of them 
in the case where a ^ —n. He considered 



^P = E' 



^A > 



A=l 



and showed that it converges for all p > 1 but that it is divergent for p < 1. 

He also gave explicit formulas for 5*2, ... 5*6 and a method (far from explicit) for 
expressing 5*^+1 as a linear combination of 5*2, ... 5*^-1. In ([3^ ) below we give a 
convolution formula for this task. 



The function w = i-Fi(a; b; z) satisfies 

zw" + {h — z)w' — aw = 
so u = w' /w satisfies the Riccati equation 

zu' + (6 — z)ti — a + zu^ = 0. 
From the Weierstrass product representation theorem, we get 



w = e 



az/b 



k=l 



n 1-^ 



Zk 



3^/Zfc 



(27) 
(28) 

(29) 



Differentiating 



logarithmically 0, 
w'{z) a 



u{z) 



w{z) 



CO 

fc=l 


l/Zk 
1 - Z/Zk 


1" 

Zk_ 


>;- 
fc=i 


1 r 

2fc i 


l Zk\ 


-1 



— 2^ Sk+iz , 



(30) 



fc=i 



where the interchange of orders of summation here is justified by the absolute conver- 
gence of the iterated series. From this we have 



zu'{z) = -J2 kSk+iz'', 

k=l 



and 



oo / k—1 



[uiz)]' = {a/hf - 2{a/h) J2 Sk+iz"" + E E S^+^S^-m+i 



(31) 



(32) 



fc=i 



k=2 \m=l 



Thus the equation ( p8|) becomes 

oo 

-Y.{h + k)Sk+iz'' + 

fc=l 



1 - 



2a 



E ^k+lZ 
k=l 
oo / k 



k+1 



a a 



+ /_^ I 2^ SmSk-m+2 ] Z' 
k=2 \m=2 



fc+1 



(33) 



Comparing the coefficients of z'^, k = 1,2, . . . in ( |55D we get: 

a{a — h) 



S, 



S. 



Sk+i 



62(6+1)' 

a{a — b){b — 2a) 
b^{b+l){b + 2)' 

k-l 



{b — 2a) Sk + 6 2^ SmSk-m+i 



A; = 3,4,.... (34) 



Kk + b) L 
This leads, in particular, to: 

_a{a-b)[a{a-b){5b + 6) + b^{b+l)] 
^~ 64(5 + 1)2(5 + 2)(6 + 3) ' 

etc., agreeing with the results found by Buchholz [ |10|1 . 

References 

[1] Watson G N 1944 A Treatise on the Theory of Bessel Functions, 2nd ed. (Cam- 
bridge: Cambridge University Press) 

[2] Ismail M E H and Muldoon M E 1995 Meth. Appl. Anal. 2 1 

[3] Obi E C 1980 J. Math. Anal. Appl. 77 465 

[4] Kishore N 1963 Proc. Amer. Math. Soc. 14 527 

[5] Muldoon M E and Raza A 1998 J. Phys. A: Math. Gen. 31 9327 

[6] Mercer A McD 1992 Internal J. Math, and Math. Sci. 15 319 

[7] Erdelyi A et al 1954 Higher Transcendental Functions, vol 2 (New York: McGraw- 
Hill) 

[8] Titchmarsh E C 1939 The Theory of Functions (Oxford: Oxford University Press) 

[9] Knopp K 1956 The Theory of Functions, Part II (New York: Dover) 

[10] Buchholz H 1951 Z. Angew. Math. Mech. 31 149 



